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A simple model to encapsulate the essential growth properties of the web of human sexual contacts 
is presented. In the model only heterosexual connection is considered and represented by a random 
growing bipartite graph where both male-female contact networks grow simultaneously. The time 
evolution of the model is analysed by a rate equation approach leading to confirm that male and 
female sexual contact distributions decay as power laws with exponents depending on influx and 
charisma of the sexes. 

PACS numbers: 02.50.cw, 05.40.-a, 89.75Hc. 

I. INTRODUCTION II. THE MODEL 



Yet another human interaction network! The web of 
human sexual contacts is introduced in Q , where the au- 
thors make extensive analysis of 1996 Swedish survey of 
sexual behaviour and show that the cumulative distribu- 
tions of number of sexual partners for males and females 
have power-law forms with exponents a/ ~ 2.54 for fe- 
males and am ~ 2.31 for males, during a single year and 



"/tot ~ 2.1, 



1.6 for entire life time. 



Human sexual behaviour is a subject area ^ in it- 
self and beyond the scope of this work. Here we will 
concentrate on capturing the essential details needed to 
form such a social network for which we will make use of 
the developments from the previous network models |^-||] 
and general consensus about human sexual behaviour. 

The model introduced in |^ has identified two basic 
principles that seem to govern the growth dynamics of 
human interaction networks, namely preferential attach- 
ment and continuous growth. This differs from the tra- 
ditional complex network consideration of connect- 
ing together a fixed number of network elements (nodes) 
which results in Poisson distributions of the connectivity. 

The usual approach in most complex network mod- 
els [^,^,^,|2| that have human attributes is to employ 
the preferential attachment process; an incoming node 
is more likely to connect to a node with higher connec- 
tions. This process yields power-law distributed connec- 
tivity in the network. A similar process takes place in a 
sociophysics model p^ , where probability of winning a 
fight depends on the number of victories, consequence of 
which is a self-organised hierarchy in the system. Again 
a well studied biological evolution model ||lj|, in which 
the power-law behaviour is due to increasing fitness of 
the species. 

Contrary to previous human interaction network mod- 
els, the model we consider here has two different types 
of nodes and the linking process is discriminative as well 
as preferential. In the proceeding section, the model is 
built by connecting new nodes to nodes with higher con- 
nections and analysed by a rate equation approach. The 
results and the discussions are presented in the final sec- 
tion. 



The growth process of the sexual contact network de- 
scribed by a random growing bipartite graph (illustrated 
in Fig. |l|) considers the following: 

(i) With probability q a new (no previous connection) 
female is connected to an existing male, which con- 
tributes to the number of females with contact de- 
gree I — 1 in the network and increases the degree 
k of the connecting male. 

(ii) With probability p a new male is connected to an 
existing female, which has the analogous effect as 
in (i). 

(iii) With probability r a link is created between an ex- 
isting male and an existing female, which increases 
the degrees of both sexes. 

In addition to above rules, initially there must be at least 
one connected male-female pair and only heterosexual 
contact is allowed among the sexes. Also the process 
r only counts once between the same pair; that is once 
a link between two individuals is created it remains and 
further contact between them does not increase their con- 
tact degrees. 
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FIG. 1. Schematic illustration of the growth processes of 
male-female sexual contact network. Male nodes (squares) 
labelled A to F and female nodes (circles) labelled 1 to 5. 
The numbers above the nodes are indicating the node degree 
/ for females and k for males, i.e. female node 1 has degree 2 
and male node D has degree 3. 



We define male sexual degree distribution Ak as the 
number of males with k female connections. The rate at 
which a male attracts females is proportional to his num- 
ber of connections and it grows linearly with increasing 
fc; in other words the connection probability of a male 
increases as he becomes more experienced. And simi- 
larly does the female sexual degree distribution Bi with 
attraction rate proportional to I. However, it seems nec- 
essary to consider some characteristics unique to both 
sexes so that we may underpin the social expectations 
about their sexual contacts. Here we introduce A and fi 
as the charisma parameters of males and females respec- 
tively. These parameters should not be confused with 
charm or beauty, which could be used as attractiveness 
of individuals. Here we are considering only the general 
attractiveness of the sexes, which are shaped by social 
expectations. 

The charisma parameters are added into the rate of 
connections given fc -I- A for males and I + fj, for females. 
Since we need fc + A > and Z + /i > for the network 
to grow, A and n are considered to be positive integers 
and small compared to the contact degrees of fc and I. 
Also, the general consensus is that males tend to exag- 
gerate their number of sexual partners. This desire puts 
females in demand and increases female charisma. Con- 
versely, the social expectation of females restricts them 
to report less and partly because of this they are more 
selective about their sexual partners. Using this infor- 
mation we argue that the rate at which females attract 
sexual partners should be greater than that for males and 
therefore fi > X. 

In this model we are considering a large system where 
the random fluctuations in the linking process can be ne- 
glected when compared with the collective behaviour of 
the individuals. Therefore it is appropriate to use mean- 
field approximation. Then we can write the time evo- 
lution of the average number of males with fc contacts 
Ak{t) as 



dAkjt) 
dt 



Ma 



[{k-l + X)Ak-i{t)-ik + X)Akit)] 



(1) 



and similarly for females with / contacts, 
dBi{t) p + r 



dt 



Mb 
f<?<5/i 



[{l^l + fi)Bi^i{t)~{l + fi)Bi{t)] 



(2) 



Since the structure of Eq. (|l|) and Eq. 1^ are the same we 
will only analyse the former equation and apply results 
to the latter. 

The first term in the square brackets indicates the con- 
tribution to Ak when a male of fc — 1 connections acquires 
a new contact, the corresponding losses are given by the 
second term. The growth of the network due to influx is 
expressed by the last term, with p being a probability of 



addition of males with fc = 1 connection. The multiplica- 
tive factor Ma is to ensure the appropriate normalization 
and is given by 



Ma{t) = J2(k + X)Ak{t). 



(3) 



The total number of males in the system is A{t) — 
J2k^k{t) and the total number of females B{t) ~ 
J2i Bi{t). From initial conditions we have A{t) — 1 + pt 
and similarly B{t) = I + qt making an assumption that 
there is one of each sex representatives before the growth 
commences. In particular, the bipartite structure ^5 16 
of the graph ensures that the total number of connections 
in the system is the same for both males and females, 
which means that L{t) = K{t) ~ J2k ^^k{t) = t + 1, 
with one link before the first time step. Also, the net- 
work does not contain monogamous individuals after the 
first time step. 

We start our analysis by looking at the moments of 
Ak{t), which are defined as 



M,{t)=Y,k'Ak{t). 



For the first few moments it is easy to show that 



and 



Mo{t) = Ma{0) + {q + r + p)t 



Mi{t) ^ Mi{0) + pt. 



(4) 



(5) 



(6) 



For large times the initial values of the moments become 
irrelevant so that we get 



Ma{t) = [q + r + p{l + X)]t, 



(7) 



which is a linear function of time and indicates that 
Eq. (|l|) is also linear in time. Using this input we can 
write Ma(t) = mat implying rria ^ q + r + p{l + A) and 
Ak{t) = akt. Substituting these relations in Eq. (|l]) leads 
to a time independent recurrence relation 

+ fc + a) flfc - (fc - 1 + X)ak-i + ^pSki- (8) 
\q + r J q + r 



Iterating Eq. (|^) we obtain 

r(fc + A) Tj^ + X + l) 



flfe 



r(^ + fc + A+l) r(A) q + r 



p. (9) 



Using the former relation for the case of females, by sim- 
ply replacing the relevant parameters gives 

rq-f/i) + + 

'' = r(^ + ^ + M + i) m JT-r'- ^''^ 

For large fc and I Eq. (||) and Eq. ( [lO[ ) are asymptotically 
equivalent to 



2 



flfe 



and 



r(i+^). 



(11) 



(12) 



The distributions ~ k^'^'^ and 6; ~ Z"'*'^ scale as 
power-laws, which is a vulnerable structure for social 
networks as far as the epidemics and the spreading of 
computer viruses are concerned [p7[-p^ . The scaling ex- 
ponent 



P(l + A) 

l~p 



depends upon male charisma parameter A and the rate 
p of addition of males into the network. Similarly, the 
exponent 



7/ = 2- 



g(l + M) 

i-g 



(14) 



depends upon female charisma parameter /i as well as 
the rate of arrival of females. However, since we have the 
constraint p + q + r = 1 implying one link at each time 
step, the growth of one network is effected by the growth 
of the other. 



III. DISCUSSION AND CONCLUSIONS 

From the definition of the model, we have four free 
parameters, 0<p<l, 0<g<l, /i>A>0 and the 
fifth one is given by the relation (p + q + r) = 1. 

Starting with a simple example, let us ignore the pro- 
cess r and take p + q — 1. This yields 



If 



g(l + M) 

1-g 



and 



Ir, 



= 2 



(l-g)(l + A) 



(15) 



We get 7/ — > 2, 7,„ — ^ oo as g ^ and 7/ — > 00, 7„j 2 
as g ^ 1. Although this is an unrealistic assumption, 
it simply illustrates the strong coupling of the influxes 
of males and females into the system. A similar strong 
coupling is also observed in where a more general 
coupled growing network is considered. 

If we pick A = 2, /i = 3 and substitute into the scaling 
exponents in Eq. (^3|) and Eq. (|l^) respectively we get 
p « 0.30 and o « 0.28 for the corresponding power-law 
exponents in ||l|]. Where 7™ = l+Om = 3.31 and 7/ — 1+ 
ctf — 3.54 for the short time behaviour analysis. In this 
example choice of values for A and /i allowed us to obtain 
p K q meaning that the increase of links due to new males 
and females is approximately the same and about 40 % of 
the links are created among the pre-existing males and 
females. The average number of partners of the sexes 
fc^fc/ ^fe = 1/p ~ 3 is the same since we have a 
close system. 

Let us now consider the distribution of contacts dur- 
ing the life time of the sexes. Here, we introduce another 



parameter /3 to the attachment rates to describe bad ex- 
periences, aging or inactivity of the sexes and hence the 
new parameter is a negative input. It can be in the range 
[0, k + \) or [0, ^ -I- ^), which seems reasonable to assume 
for k,l > 2. We can denote Ai = A — (3m and jii — ji — (3f 
where we assign a unique /3 parameter to each kind to 
discriminate their sexual behaviour. The new scaling ex- 
ponents 



(13) and similarly 



If 



7n 



g(l + A^i) 



P(l + Ai) 

1 - p ' 



(16) 



(17) 



Using the same values we used earlier for the free param- 
eters gives 7/ = 1 + a/4„t — 3.1 for fif ^ 1.17 the same 
for all females and 7™ = 1 + otmtot — 2-6 for ~ 1-6 for 
all males. Consequently we have /?/ < /3m implying that 
the sexual connections of males tend to suffer more from 
aging, bad experiences or inactivity than that of females. 
Perhaps this may be the reason why males resort to lies 
to maintain self image. 

In conclusion, our simple model seems adequate to sim- 
ulate the dynamics of the web of human sexual contacts. 
The gratifying feature of rate equation approach is that it 
gives flexibility to include diverse range of parameters to 
model any situation within the frame of this work. The 
charisma parameters used here are shaped by the western 
social expectations and they can be changed to reflect the 
views of other societies. Also, one can introduce fitness, 
multiplicative quenched disorder into the system as in 
|]l2| to account the activity level of different age groups. 
Finally, a challenging issue to address would be to study 
male-female degree correlations of this model. 
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